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Abstract 

Using 5.8x 10^ J/V' events collected in the BESII detector, the radiative decay J/if) ^ 
700 — >■ ^K'^K'K^K^ is studied. The cfxp invariant mass distribution exhibits a 
near-threshold enhancement that peaks around 2.24 GeV/c^. A partial wave anal- 
ysis shows that the structure is dominated by a ^ state (?7(2225)) with a mass of 
2.24l[j:[J^l[j:[|^ GeV/c2 and a width of 0.19±0.03lg;[j^ GeV/c^. The product branch- 
ing fraction is: Br(J/^ ^ 7r/(2225)) ■ Br(r/(2225) ^ #) = (4.4± 0.4± 0.8) x lO"''. 



1 Introduction 



In addition to the spectrum of ordinary light qq meson states, QCD-motivated 
models predict a rich spectrum of gg glueballs, qqg hybrids and qqqq four quark 
states [1]. Radiative J/'ip decays provide an excellent laboratory for testing 
these predictions, and systems of two vector particles have been intensively 
examined for signatures of gluonic bound states [2]. Pseudoscalar enhance- 
ments in pp and uouo final states have been seen in radiative J /ip decays [3-7]. 
Recently, a near-threshold scalar, the /o(1810) or X(1810) was reported in 
the oj(f) invariant mass distribution from the doubly OZI suppressed decays of 
J/%1) ^ ju!(l) [8], thereby adding an additional puzzle to the already confusing 
spectrum of low-lying scalar mesons [9,10]. 

Structures in the (jxp invariant-mass spectrum have been observed by several 
experiments both in the reaction 7r"p — > (fxpn [11] and in radiative J/ip de- 
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cays [12-14]. The ?7(2225) was first observed by the MARK III collaboration in 
J /ijj radiative decays J /ip ^ 700 [12]. A fit to the 00 invariant-mass spectrum 
gave a mass of 2230 ± 25 ± 15 MeV/c^ and a width of 150t^[]° ± 60 MeV/c^ 
The production branching fractions are Br{J/'4) 777(2225)) • Br{r]{2225) 
(t)(t)) = (3.3 ± 0.8 ± 0.5) X 10-^ for the -fK^R-K+R- mode and Sr( J/'0 ^ 
777(2225)) •Sr(r7(2225) ^ 00) = (2.7±0.6±0.6) x 10"^ for the 7K+/^-/r^ir£ 
mode. An angular analysis of the structure found it to be consistent with a 
^ assignment. It was subsequently observed by the DM2 collaboration, also 
in J/'^ — > 700 decays [13,14]. 

In this letter we present results from a high statistics study of J/%1) ^ 700 in 
the ^R^R~RgR^ final state, using a sample of 5.8 x 1^'^ J/ip events collected 
with the BESII detector at the Beijing Electron- Positron Collider (BEPC). 
The presence of a signal around 2.24 GeV/c^ and its pseudoscalar character 
are confirmed, and the mass, width, and branching fraction are determined by 
a partial wave analysis (PWA). 



2 BES detector and Monte Cetrlo simulation 



BESII is a large solid-angle magnetic spectrometer that is described in detail 
in Ref. [15]. Charged particle momenta are determined with a resolution of 
Gp/p — 1.78%-\/l (with p in GeV/c) in a 40-layer cyhndrical main drift 
chamber (MDC). Particle identification is accomplished using specific ioniza- 
tion (dE/dx) measurements in the MDC and time-of-flight (TOF) measure- 
ments in a barrel-like array of 48 scintillation counters. The dE/dx resolution 
is <y dE/dx — 8.0%; the TOF resolution is gtof = 180 ps for the Bhabha events. 
Outside of the TOF counters is a 12-radiation-length barrel shower counter 
(BSC) comprised of gas tubes interleaved with lead sheets. The BSC measures 
the energies and directions of photons with resolutions of (Je/E ~ 21%/V^ 
(with E in GeV), cr^ = 7.9 mrad, and = 2.3 cm. The iron fiux return of the 
magnet is instrumented with three double layers of counters that are used to 
identify muons. 

In this analysis, a GEANT3-based Monte Carlo (MC) simulation program 
(SIMBES) [16], which includes detailed consideration of the actual detector 
responses (such as dead electronic channels), is used. The consistency between 
data and Monte Carlo has been checked with many high-purity physics chan- 
nels [16]. 
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3 Event selection 



Since the is difficult to identify in BESII, its detection is not required in 
the selection of J/ if) -fK+R-R^sKl events. J/i/j -fK+K-n+n-{Kl)miss 
candidates arc selected from events with four charged tracks with net charge 
zero in the MDC and with one or two isolated photons in the BSC. Charged 
tracks are required to be well fitted to a helix, be within the polar angle region 
I cos ^1 < 0.8, and have atransverse momentum larger than 50 MeV/c. For each 
track, the time-of-flight (TOF) and specific ionization (dE/dx) measurements 
in the MDC are combined to form a particle identification confidence level 
for the TT, K and p hypotheses; the particle type with the highest confidence 
level is assigned to each track. The four selected charged tracks are required to 
consist of an unambiguously identified K^, K~ , tt"*" and 7r~ combination. Each 
candidate photon is required to have an energy deposit in the BSC greater 
than 60 MeV, to be isolated from charged tracks by more than 20°, have an 
angle between the cluster development direction in the BSC and the photon 
emission direction less than 30°, and have its first hit in the BSC within the 
first six radiation lengths. A one-constraint (IC) kinematic fit is performed with 
the J/V' — > ^K^K~7r'^7r'{Kl)^^^^ hypothesis. For events with two isolated 
photons, the one that gives the smallest x^{lK~^K~7[~^n~{Kl)^^g^) is selected. 
In order to improve the K^K^ mass resolution, a 2C-kinematic fit is performed 
by adding a Kg mass constraint. Events with Xicil^'^^'^^^'i^Dmiss) 
20 and xlch^+K- K^si^D^iJ < 20 are retained. 

The KgKl versus K'^K~ invariant masses are plotted in Fig. 1 (a), where 
a cluster of events corresponding to 00 production is evident. Because the 
processes J/ip (fxf) and J/ih -R^cfx}) are forbidden by C-invariance, the 
presence of two 0's is a clear signal for the radiative decay J/'0 — >• 700. 
The histogram in Fig. 1 (b) shows the K^K" invariant mass distribution 
after the requirement that the K^K^ invariant mass is inside the — > K^K^ 
signal region {\M{KgK^) —m^\ < 0.025 GeV/c^). The histogram in Fig. 1 (c) 
shows the KgKl invariant mass distribution after the requirement that the 
K^K~ invariant mass is inside the — > K^K^ signal region {\M[K^ K~) — 
< 0.0125 GeV/c^). The histogram in Fig. 2 (a) shows the K+K-K%Kl 
invariant mass distribution for events where the K^K~ and K^K^ invariant 
masses lie within the 00 mass region {\M{K^ K^) — m^\ < 0.0125 GeV/c^ 
and \M{K^sKl) - < 0.025 GeV/c^). There are a total of 508 events, 
which survive the above- listed criteria (optimized for low 00 masses), with 
a prominent structure around 2.24 GeV/c^. The phase space invariant mass 
distribution and the acceptance versus 00 invariant mass are also shown in 
Fig. 2 (a) as the dashed histogram and dotted curve, respectively. The peak is 
also evident as a diagonal band along the upper right-hand edge of the Dalitz 
plot, shown in Fig. 2 (b). The asymmetry in the dalitz plot of data is caused 
by detection efficiency. 
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Fig. 1. (a) The K^K^ versus K 

invariant mass distribution when the K^K^ invariant mass is in the 



signal region. The arrows show the (j) K^K^ signal region, (c) The K^K^ invari- 
ant mass distribution when the K^K" invariant mass is in the (j) K^K~ signal 
region. The arrows show the (j) — > K^K^ signal region. 



Non-(^(/) backgrounds are studied using events in the sideband regions shown 
in Fig. 1 (a). Figure 3 (b) shows the K~^K~KgKl invariant mass of events 
within the K+K' sideband region (1.069 GeV/c^ < M{K+K-) < 1.094 
GeV/c2 and \M{K'^sK'l) - < 0.025 GeV/c^), and Fig. 3 (c) shows the 
corresponding spectrum of events within the —>■ K^K^ sideband region 
{\M{K+K-)-m^\ < 0.0125 GeV/c^ and 1.082 GeV/c2< M{K'^sK) < 1-132 
GeV/c^). Figure 3 (d) shows the events in the corner region, which is de- 
fined as (1.069 GeV/c2< M{K+K-) < 1.094 GeV/c^ and 1.082 GeV/c2< 
M{KgKi^) < 1.132 GeV/c^). The background, estimated by summing up the 



5 




5" 5 



o > 
c: CD 



LU 



T 



3 - 



2.25 2.5 2.75 



1 ^ 
1 



(a) (GeV/c'' 



(b)M'('H)) (^^KgK^) (GeVV) 



Fig. 2. (a) The K'^K^K'^K^ invariant mass distribution for J/ip ^ ^cpcp candidate 
events. The dashed histogram is the phase space invariant mass distribution, and 
the dotted curve indicates how the acceptance varies with the (f)(j) invariant mass, 
(b) Dahtz plot for J/ijj —> jcjxj) candidate events. 



scaled event yields in Figs. 3 (b) and (c) and subtracting that in Fig. 3 (d), is 
shown as the dashed histogram in Fig. 3 (a). No sign of an enhancement near 
the mass threshold is evident in the non-(f)(f) background events. 
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Fig. 3. The K'^K~KgKl invariant mass distributions for (a) events in the (fx/) 
signal region; (b) events in the — K'^K" sideband region; (c) events in the 
(j) KgK^ sideband region; (d) events in the corner region, as described in the 
text. The dashed histogram in (a) shows the background distribution obtained from 
the sideband evaluation. The dotted histogram in (a) shows the (jxf) invariant mass 
distribution of inclusive J/i/^ Monte Carlo samples 
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The backgrounds in the selected event sample are studied with Monte Carlo 
simulations of a number of potential background decay channels listed in the 
PDG tables [17]. The main background originates from J/'0 — > (j)K*K, K* 
Kn^ (+ c.c). Using the PDG's world average branching fraction for this mode, 
we estimate that about 40 events from this channel are in the (fxp invariant 
mass signal region. However, the simulation also shows that they do not peak 
at low (pcf) masses. Using a Monte Carlo sample of 35. IM inclusive J/ip decay 
events generated with the LUND-charm model [18], we searched for other 
possible background channels. None of the simulated channels produce a peak 
at low 00 invariant masses. The dotted histogram in Fig. 3 (a) shows the 00 
invariant mass distribution and the normalized background estimated with 
inclusive J/ip Monte Carlo samples where the events oi J/ijj 700 and 
J/ip —>■ ■yr]c{r]c 4>4>) are removed. 



4 Partial wave analysis 

A partial wave analysis (PWA) of the events with M(00) < 2.7 GeV/c^ was 
performed. The two-body decay amplitudes in the sequential decay process 
J/ip ^ jXjX — i> 00, — > K^K~ and K^K^ are constructed using the 
covariant helicity coupling amplitude method [19]. The intermediate resonance 
X is described with the normal Breit-Wigner propagator BW = 1/(M^ — 
s — IMF), where s is the 00 invariant mass-squared and M and F are the 
resonance's mass and width. The amplitude for the sequential decay process 
is the product of all decay amplitudes and the Breit-Wigner propagator. The 
total differential cross section da/d^ is 



where A{J^'~') is the total amplitude for all resonances whose spin-parity are 
J^*", and BG denotes the background contribution, which is described by a 
non-interfering phase space term. 

The relative magnitudes and phases of the amplitudes are determined by an 
unbinned maximum likelihood fit. The basis of likelihood fitting is the calcu- 
lation of the probability that a hypothesized probability distribution function 
would produce the data set under consideration. The probability to observe 
the event characterized by the measurement ^ is 



da 



|A(0-+) + A(0++) + ^(2++)|2 + BG, 



(1) 



P{i) 



(2) 
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where = and e(^) is the detection efficiency. The normahzation integral 
/ d^uj{^)e{^) is done with a weighted phase space MC sample; the details are 
described in Ref. [20]. The joint probability density for observing the events 
in the data sample is 



N N 



^^n-(«=njf|i5. (3) 

For technical reasons, rather than maximizing jC, S — —InjC is minimized, i. e. 



For a given data set, the second term is a constant and has no impact on the 
determination of the parameters of the amplitudes or on the relative changes 
of <S values. So, for the fitting, ln>C, defined as: 



is used. The free parameters are optimized by MINUIT [21]. In the minimiza- 
tion procedure, a change in log likelihood of 0.5 represents a one standard 
deviation effect for the one parameter case. 

For the production of a pseudoscalar, only V waves are allowed in both the 
radiative decay J/ip ^ 7X and the hadronic decay X (jxf). For the case of 
a scalar, both S and V waves are possible in both the radiative and hadronic 
decays; only the S wave is considered in the fit. For the production of a 2+ 
resonance, there are five possible amplitudes for both the radiative and the 
hadronic decays, one S wave, three V waves and one (L = 4) wave. In this 
case, only S and V waves in both decays, corresponding to the lower overall 
spin of the 00 system, are considered. (The (L = 4) wave is ignored in this 
analysis.) 

When J/'^ — > jX, X — > 00 is fitted with both the 00 and 7X systems 
in a P wave, which corresponds to a X = ^ pseudoscalar state, the fit 
gives 195.5111? events with mass M = 2.24l[J;gi GeV/c^, width F = 0.19 ± 
0.03 GeV/c^, and a statistical significance larger than 10 a. The errors are 
statistical only. Using a selection efficiency of 3.29%, which is determined from 
the Monte-Carlo simulation using the magnitudes and phases of the partial 
amplitudes from the PWA, we obtain a product branching fraction of: 

Br{J/ij 77^(2225)) ■ Er(r/(2225) ^ 00) = (4.4 ± 0.4) x 10'^. 
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Details of the fitting procedure and the detection efficiency determination 
can be found in Ref. [20]. Figure 4 (a) shows a comparison of the data and 
MC projections of the 00 invariant mass distribution for the fitted parameters. 
Comparisons of the projected data and MC angular distributions for the events 
with (f)(f) invariant mass less than 2.7 GeV/c^ are shown in Figs. 4 (b)-(f). 




Fig. 4. Comparisons between data and MC projections for the (fxl) invariant mass 
distribution and the angular distributions for events with 00 invariant mass less 
than 2.7 GeV/c^ using the fitted ?7(2225) resonance parameters. The points with 
error bars are data, the solid histograms are the MC projections, and the dashed 
lines are the background contributions, (a) The 00 invariant mass distribution; (b) 
the polar angle of the radiative photon (O^f); (c) the polar angle of the in the 00 
rest system {6^); (d) the polar angle of the kaon in the rest system {Or)', (e) the 
azimuthal angle of the in the 00 rest system (00); and (f) the x distribution, 
where % is the azimuthal angle between the normal directions of the two decay 
planes in the 00 rest frame. 



We also tried to fit the resonance with 0"^+ and 2++ spin-parity hypotheses 
using all possible combinations of orbital angular momenta in the 00 and 

7(00) systems. The log-likelihood values of the best fits arc worse than that 
of the 0^+ assignment by 95 and 27 for the O"*"^ and 2++ assignments, respec- 
tively. We therefore conclude that the J^*" of the resonance strongly favors 

0-+ (77(2225)). 

Because of the possible existence of a coherent J/ip ^ '-/(pel) phase space am- 
plitude, we also fitted with an interfering phase space (0^) term included. The 
log-likelihood value S improves by 0.4, which suggests that the contribution 
from 0~+ phase space is neghgible. 
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If an additional resonance is included in the fit, the significance of the addi- 
tional resonance after reoptimization is 0.8 a. 2.1 a and 3.3 a for the 0^+, O^"*" 
and 2+"'" assignments of the additional resonance, respectively. The differences 
between the results including and not including the additional 0+"'" or 2++ are 
included in the systematic errors; the systematic uncertainty contributions for 
the mass, width, and branching fraction of the 77(2225) are +0.01 GeV/c^, 
+0.04 GeV/c^ and +2.6%, respectively. 



5 Systematic error 

The systematic uncertainties are estimated by considering the following: the 
uncertainties in the modeling of the background, different Breit-Wigner pa- 
rameterizations, the possible presence of additional resonances, the simulation 
of the MDC wire resolution, as well as possible biases in the fitting procedure. 
The uncertainties in the background include the uncertainty in the treatment 
of the background in the fitting. We also tried to subtract the backgrounds 
determined from the sidebands in the fit, and the differences are taken as sys- 
tematic errors. Since the enhancement is in the near-threshold region of the 
invariant mass spectrum, a fit using a Breit-Wigner with a momentum 
dependent width [22] is also performed, and the differences between the fit 
with a constant width Breit-Wigner are included as systematic errors. Possi- 
ble fitting biases arc estimated from the differences obtained between input 
and output masses and widths from Monte Carlo samples, which are generated 
as J/ip ^ 70""'', 0""*" — >• 4>(j) using the fitted parameters. The total systematic 
errors are obtained by adding the individual errors in quadrature. The total 
systematic errors on the mass and width are determined to be lo!o2 GeV/c^ 
and lo!o4 GeV/c^, respectively. 

For the systematic error on the branching fraction measurement, the system- 
atic uncertainties of the photon detection efficiency and the particle identifi- 
cation efficiency, as well as the and Kg decay branching fractions, the mass 
and width uncertainties of ?7(2225), and the total number of J/ip events [23] 
are also included. The total relative systematic error on the product branching 
fraction is 



6 Summciry 

Using 5.8 X lO'^J/ip events measured in the BESll detector, we studied the 
radiative decay J/ip ^ ^(jxj) ■jK'^K^KgKl. A structure (?7(2225)) in the 
near-threshold region of (fxf) invariant mass spectrum is observed. 
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A PWA shows that the structure is dominated by a ^ state with a mass 
2.2Ato^ito'^i GeV/c^ and a width O.lQiO.OStn! GeV/c\ The product branch- 
ing fraction is measured to be: 

Br(J/jjj 777(2225)) • Br(ri(2225) 00) = (4.4 ± 0.4 ± 0.8) x 10"^ 
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